Adiabatic processes need not correspond to optimal work 
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The minimum work principle states that work done on a thermally isolated equilibrium system 
is minimal for the adiabatically slow (reversible) realization of a given process. This principle, one 
of the formulations of the second law, is studied here for finite (possibly large) quantum systems 
interacting with macroscopic sources of work. It is shown to be valid as long as the adiabatic energy 
levels do not cross. If level crossing does occur, counter examples are discussed, showing that the 
minimum work principle can be violated and that optimal processes are neither adiabatically slow 
nor reversible. 



The second law of thermodynamics^*^, formulated 
nearly one and half century ago, continues to be un- 
der scrutinity^i^. While its status within equilibrium 
thermodynamics and statistical physics is by now wcU- 
settlediiSi^, its fate in various border situations is far from 
being clear. In the macroscopic realm the second law is 
a set of equivalent statements concerning quantities such 
as entropy, heat, work, etc. However, in more general 
situations these statements need not be equivalent and 
some, e.g. those involving entropy, may have only a lim- 
ited applicabilitji^. In contrast to entropy, the concept 
of work has a well-defined operational meaning for finite 
systems interacting with macroscopic work sourceSiS,. It 
is, perhaps, not accidental that Thomson's formulation 
of the second law^'^'^ — no work can be extracted from 
an equilibrium system by means of a cyclis process — 
was proveni^ both in quantum and classical situation. 

Here we study the minimum work principle which 
extends the Thomson's formulation to non-cyclic pro- 
cesseai'^'^, and provides a recipe for reducing energy 
costs. After formulating the principle and discussing it 
for macroscopic systems, we investigate it for finite sys- 
tems coupled to macroscopic sources of work. Its domain 
of validity there is found to be large but definitely lim- 
ited. These limits are illustrated via counterexamples. 

The setup. Consider a quantum system S which is 
thermally isolatedliSi^: it moves according to its own dy- 
namics and interacts with an external macroscopic work 
source. This interaction is realized via time-dependence 
of some parameters R{t) = {i?i(t), i?2(i), ...} of the sys- 
tem's Hamiltonian H{t) = H{R{t)}. They move along 
a certain trajectory R(t) which at some initial time U 
starts from Ri — R{ti), and ends at i?f = R{ti). The ini- 
tial and final values of the Hamiltonian are Hi = H{Ri} 
and Hf — H{Rf}, respectively. Initially S is assumed to 
be in equilibrium at temperature T = 1//3 > 0, that is, 
S is described by a gibbsian density operator: 

= exp(-/3i7i)/^i, Zi = tTe-^"\ (1) 

As usual, this equilibrium state is prepared by a weak 
interaction between S and a macroscopic thermal bath 
at temperature Ti2i&, and then decoupling S from the 
bath in order to achieve a thermally isolated processiiSi^. 



The Hamiltonian H{t) generates a unitary evolution: 

th^^pit) = [H{t),p{t)], p{t) = U{t)p{ti)U\t), (2) 

with time-ordered U{t) — cxp[—\ dsH{s)]. The work 
W done on S readsi*^*^ 

j^ Attv [p[t)H{t)] = ir[Hip[U)] - tr[Fip(ii)], (3) 

where we performed partial integration and inserted 1^)1. 
This is the energy increase of S, which coincides with the 
energy decrease of the source. 

The principle. Let S start in the state and let R 
move between R{ and i?f along a trajectory R{t). The 
work done on S is W . Consider the adiabatically slow 
realization of this process: R proceeds between the same 
values Ri and i?f and along the same trajectory, but now 
with a homogeneously vanishing velocity, thereby taking 
a very long time tf — ti, at the cost of an amount work 
W . The minimum- work principle then assertsisSi^ 

W>W. (4) 

This is a statement on optimality: if work has to be ex- 
tracted from S, W is negative, and to make it as negative 
as possible one proceeds with very slow velocity. If during 
some operation work has to be added {W > 0) to S, one 
wishes to minimize its amount, and operates slowly. For 
thermally isolated systems, adiabatically slow processes 
are reversible. This is standard if S is macroscopioiiSi^, 
and below it is shown to hold for a finite S as well, where 
the definition of reversibility extends unambiguously (i.e., 
without invoking entropy)'^. 

In macroscopic thermodynamics the minimum work 
principle is derived^'^ from certain axioms which ensure 
that, within the domain of their applicability, this princi- 
ple is equivalent to other formulations of the second law. 
Derivations in the context of statistical thermodynamics 
are presented inSiSiili. We discuss one of them now. 

The minimal work principle for macroscopic systems 
is proven in two steps: first one considers the relative 
entropy tr[p(tf ) In p(tt) — p{tf) In peq{-fff)] between the fi- 
nal state p(if) given by ^ and an equilibrium state 
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Poq{Hf) = exp{—PH{)/Zf, Z{ = tve~^^^, a state corre- 
sponding to the final Hamiltonian Hi and the same tem- 
perature T ^ 1/(3. As follows from tr[p(if) In p{t{)] = 
tr[/9(ii) ln(o(ti)]. This combined with and the non- 

negativity of relative entropy^ yields: 

W > F{Hf) - F{Hi) = rintre"'^^' - Tlntre"'^^', (5) 

where F{Hi) and F{Hf) are the gibbsian free energies 
corresponding to p{ti) and pcq{H{), respectively. 

There are several classes of macroscopic systems for 
which one can show that the free energy difference in the 
RHS of (0 indeed coincides with the adiabatic worliS*2iiS. 

Finite systems. For an arbitrary iV-level quantum sys- 
tem S, Eq. © does not have the needed physical mean- 
ing, since in general F{Hf) — F{Hi) does not coincide 
with the the adiabatic work. It is know»ii that for finite 
systems the final density matrix p{t{) given by (0 need 
not coincide with pcq{Hi) = exp{—(3H{)/Zi. This fact 
was recently applied for certain irreversible processesi^. 

Thus we need an independent derivation of I^J. Let 
the spectral resolution of H{t) and p{ti) be 
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At t = ti we order the spectrum as 

£l{ti) < ... < SNiti) => Pi> ... > PN- 
For ti < t < tf we expand on the complete set |n, t): 

N 

u{t)\k, ti) = w I"' 

n=l 

and use © to obtain: 
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A similar formula can be derived to express the adia- 
batic work W in coefficients afc„(tf). From the definition 
|afe„(tf)p = |(n,if|t/|fc,ti)|2 it follows that 
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With help of the identity: ^ 
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J2Z=i[^rn+i-£7n] Y.n=i ^n, wc obtaiu using JTHl Ell the 
general formula for the difference between non-adiabatic 
and adiabatic work: 
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[£m+l(^f) - em(if)] ©m, (12) 
N 

EE^^^(i«^"(^f)i'-i«'=«(*f)i')- (13) 



To understand the meaning of this formula, let us first 
assume that the ordering ||SJ) is kept at t — tf-. 



£l{tf 



< 



< £jv(tf). 



(14) 



If different energy levels did not cross each other (and 
equal ones do not become different), Eq. p4|l is im- 
plied by Eq. ©. According to non-crossing rule^'^, if 
H{R} is real and only one of its parameters is varied 
with time, H14|l is satisfied for any discrete-level quan- 
tum system: level-crossing, even if it happens in model- 
dependent calculations or due to approximate symmetry, 
does not survive arbitrary small perturbation where it is 
substituted by avoided crossing (for a more general H{R} 
the conditions prohibiting level-crossing are more restric- 
tive; seei^). No level-crossings and natural conditions of 
smoothmess of H(t) are sufficient for the standard quan- 
tum adiabatic theorem-'^'* to ensure 



afen(tf) = Skn- 

Combined with ijHlEll, Eq.(23Il brings 
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Eqs. (Uni 113 Ell) together with 6™ > extend the mini- 
mum work principle to cases where the adiabatic work 
is not equal to the difference in free energies. 

Level crossing. The above non-crossing condition 
raises the question: Is the minimum work principle also 
valid if the adiabatic energy levels cross? Before address- 
ing this question in detail, let us mention some popular 
misconceptions which surround the level-crossing prob- 
lem: 1) The no-crossing rule is said to exclude all cross- 
ings. This is incorrect as the exclusion concerns situa- 
tions where, in particular, only one independent parame- 
ter of a real Hamiltonian H{R} is variedtf.. Two param- 
eters can produce robust level-crossing for such Hamilto- 
nians. 2) It is believed that once levels can cross, Ae 0, 
the very point of the adiabatic theorem disappears as the 
internal characteristic time h/ Ae of S is infinite. This 
view misidentifies the proper internal time as seen be- 
low; see alsol^ in this context. 3) It is sometimes believed 
that crossing is automatically followed by a population 
inversion. We shall find no support for that. 

As a first example we consider a spin-1/2 particle with 
Hamiltonian 



H{t) = hi{s)ai ~ h3{s)a3, 



t/T, 



(16) 



r! = l k=l 



where ci, and (72 = ificra are Pauli matrices, and 
where s is the reduced time with t being the character- 
istic time-scale. The magnetic fields hi and /13 smoothly 
vary in time. Assume that i) for s — s- Sj < and for 
s ^ Sf > 0, hi{s) and h^ls) go to constant values 
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sufRciently fast; ii) for s ^ one has: hi{s) ~ ais^, 
/i3(s) ~ —0:35, where ai and are positive constants. 
in) hi{s) and /i3(s) are non-zero for all s, Si < s < S{, 
except s = 0. Not all these points are needed, but we 
choose them for clarity. Generalizations are indicated 
below. One writes (|16|l as 



H ^ 



-h^is) hi{s) 
hi{s) /i3(s) 



ei(s) 



cost 
sin( 



sin 9{s) 
■ cos 6(s) 



where 9{s) = — arctan[/ii(s)//i3(s)] is a parameter in the 
interval — 7r/2 < 9 < 7r/2, and where 



e,{s) = sg{s)^hlis) + hlis), e2(,s) = (17) 

are the adiabatic energy levels which cross at s = d{s) — 
(y^^ is defined to be always positive). The adiabatic 
eigenvectors are, H{s)\k,s) — ek{s)\k,s), k = 1,2, 



cos ^9{s) 
sin h9{s) 



12, s) = 



-sini6l(s) 
cos h9{s) 



.(18) 



Both the eigenvalues and the eigenvectors are smooth 
functions of s. Eq. (|14|) is not valid, and (|l()m3(l imply: 



W -W ^ -2^hl{si) + hl{si)ei, TSi^tf. (19) 

Naively this already proves the violation. More care- 
fully, our strategy is now to confirm H15(l in the adiabatic 
limit r — > cx) and thus to confirm that Qi > 0, implying 
that the minimum work principle is indeed violated. To 
this end we apply the standard adiabatic perturbation 
theory Substituting jnj into |5J) one has: 



O-kn 



N 

E 

m— 1 



afe™(t)e^'''''^*'[""^*'^~"'"^*'^l {n\dt\m). (20) 



As |1) and |2) in (|18|) are real, {n\n) = 1 implies 
{n\dt\n) = 0. Since {n\dt\m) = ^{n\ds\m) the RHS of 
(|20|l contains a small parameter 1/t. It is more conve- 
nient to introduce new variables: akn{i) = ^kn + bkn{t): 
bkn{ti) — 0. To leading order in 1/r, bkn can be neglected 
in the RHS of ((JOl, yielding for ak^n{ti) = bk^niU)- 



\ak^n{ti)Y 



,(21) 



while |afefe(tf)|2 = 1 - E„5^fc |afcn(if)p. In ^ we put 
ST = <, s't = t'. For our model (|16II18|) . du[ei{u) - 
£2('")] — 2 duei{u) has only one extremal point, at 
s = 0. We also have from (|18() 
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For large r the integral in (|21|l can be calculated with 
use of the saddle-point method: 



T 



(2|9.|1)W 



hi 



Aral ' 



(23) 



Eqs. H21I I23|l extend the adiabatic theorem H15() for the 
level-crossing situation. More general statements of sim- 
ilar adiabatic theorems can be found in Refii^. Inserting 
©1 = {Pi - P2)\ai2{U)? > in Eq. confirms the 

violation of the minimum work principle. Eq. H23|l shows 
that the role of the proper internal characteristic time is 
played by ha\/a\ rather than by h/{£i — £2). 

More generally, if \Jh\{s) + h\{s) is a smooth function 
for all real s (e.g., it is not (x |s|), there are no crossings of 
eigenvalues and (@J) is valid. If both hi{s) and ^.3(5) are 
linear for s — *■ 0, the leading term presented in H23() van- 
ishes due to (2|9s|l)^|s^o = 0, and one needs the second- 
order in the saddle-point expansion, to be compared with 
the second-order term of the adiabatic perturbation the- 
ory. This leads to the same physical conclusions as (|23|) 
did, but with |ai2(if)P t'^. 

One can calculate |afe„(tf)| yet in another limiting case, 
where the characteristic time r is very short. It is well- 
known^** that in this limit energy changes can be calcu- 
lated with frozen initial state of S. For the present sit- 
uation this leads from ifT^ to |ai2(if)P = |a2i(tf)P = 
|(l,tf I 2,ii)|2 = sm^\[9{tf) - 6'(ii)], and thus to Oi = 

(pi — P2) sin^ ^[^(^f) — ^^(^i)], again positive. 

Exactly solvable model with level crossing. Consider a 
two-level system with Hamiltonian 



H{t) = huj 



scos^s is sin 2s 



is sin 2s ssin^s 



s='-, (24) 



where r is the characteristic time-scale, and w is a con- 
stant. For Si < denote the adiabatic energy lev- 
els as £i(si) = hwsi < £2(31) = 0. They cross when 
s passes through zero. Eq. (|20|l for this model can 
be solved exactly in terms of hypergeometric functions. 
Postponing the detailed discussion, we present in Fig. 
n the behavior of |ai2(sf)p as a function of r. Since 
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FIG. 1: Amplitude |ai2(sf)p versus the time-scale t for Si = 
— 1.5, S{ = 1.5 and uj — 1. Full oscillating curve: the exact 
value which can reach unity. Dotted curve: result from a 
first-order adiabatic perturbation theory. The smooth curve 
presents a saddle-point approximation analogous to I23II . 

from 1121 124|l one has for the present model W — W = 
—hujsi |ai2(sf)p tanh(i/3?i(jj |si| ), violations of the mini- 
mum work principle exist for Sf > 0, and they are maxi- 
mal for |ai2(sf)p 1. This is seen to be the case in Fig. 
^for some r near t — 1. Notice that both the first-order 
perturbation theory and the saddle-point approximation 
are adequately reproducing |ai2(sf)p for r > 10. 
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Let S has a finite amount of levels, and two of them 
cross. For quasi-adiabatic processes (r is large but finite) 
the transition probability between non-crossing levels is 
exponentially small^^'^^, while as we saw it has power-law 
smallness for two crossing levels. Then one neglects in 
(|10|l the factors \ak^n{t{)\'^ coming from any non-crossed 
levels k and n, and the problem is reduced to the two-level 
situation. Thus already one crossing suffices to detect 
limits of the minimum work principle. The reduction to 
the two-level situation takes place also in a macroscopic 
system which has few discrete levels at the bottom of 
a continuous spectrum, since for low temperatures these 
levels can decouple from the rest of the spectrum. 

Cyclic processes and reversibility. The above results 
do not imply any violation of the second law in Thom- 
son's formulation^: no work is extracted from S during 
a cyclic process, Wc > 0. We illustrate its general proof 
in the context of the level crossing model given by 
I18|l . Assume that the trajectory R{t) — {hi{t), h2{t)) 
described there is supplemented by another trajectory 
R'{t) which brings the parameters back to their initial 
values (/ii(ti), /i3(ti) ) so that the overall process R + R' 
is cyclic. If R' crosses the levels backwards, then at the 
final time of R' Eq. is valid, and (|10H15() imply: 

Wc = \ai2\\pi - P2)[e2iU) - ei{U)] >Wc = 0, (25) 

where |ai2p < 1 now corresponds to the full cyclic pro- 
cess R + R'. Eq. H25|l confirms the intuitive expectation 
that non-adiabatic process are less optimal. In particu- 
lar, this is valid if R' is exactly the same process R moved 
backwards with the same speed. Then Wc = means 
that i? is a reversible process in the standard thermody- 
namical sensei*^*^. If R' does not induce another level 
crossing, i.e., /ii(s) and /i2(s) in Eq. p6|l return to their 
initial values without simultaneously crossing zero, then 
£i{tf) = £2{ti), £2 (if) = £i{ti) and Eqs. l(Tn|[I3) imply 

VKc = {Pi -P2)[£2{ti) " eiiU)] >Wc^ |aii|2 Wc > 0. 



In contrast to H25() , non-adiabatic processes are more op- 
timal it R + R' contains one level-crossing (or an odd 
number of them). We thus have found here a violation 
of the minimum work principle for a cyclic process. 

In conclusion, we have studied the minimum work 
principle for finite systems coupled to external sources 
of work. As compared to other formulations of the sec- 
ond law, this principle has a direct practical meaning as 
it provides a recipe for reducing energy costs of various 
processes. We gave its general proof and have shown that 
it may become limited if there are crossings of adiabatic 
energy levels: optimal processes need to be neither slow 
nor reversible. Already one crossing suffices to note vi- 
olations of the principle. If this is the case, the optimal 
process occurs at some finite speed. 

Level-crossing was observed, e.g., in molecular and 
chemical physics!^. It is not a rare effect: if the number 
of externally varied parameters is larger then two, then 
for typical spectra level crossings are even much more 
frequent than avoided crossingsii. It is possible that the 
presented limits of the minimum work principe may serve 
as a test for level crossings. 

Together with the universal validity of Thomson's for- 
mulation of the second lawSi^, the limits of the principle 
imply that the very equivalence between various formu- 
lations of the second law may be broken for a finite sys- 
tem coupled to macroscopic sources of work: different 
formulations are based on different physical mechanisms 
and have different ranges of validity. Similar results on 
non-equivalence of various formulations of the second law 
were found in Refi^, where for a quantum particle cou- 
pled to a macroscopic thermal bath, it was shown that 
some formulations, e.g., Clausius inequality and positiv- 
ity of the energy dispersion rate, are satisfied at suffi- 
ciently high temperatures of the bath, but can be invalid 
at low temperatures, that is, in the quantum regime. 
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